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Abstract 

We investigate the universal property of curvatures in surface models which 
display a flat phase and a rough phase whose criticality is described by the 
Gaussian model. Earlier we derived a relation between the Hessian of the 
free energy and the Gaussian coupling constant in the six- vertex model. Here 
we show its validity in a general setting using renormalization group argu- 
ments. The general validity of the relation is confirmed numerically in the 
RSOS model by comparing the Hessian of the free energy and the Gaussian 
coupling constant in a transfer matrix finite-size-scaling study. The Hessian 
relation gives clear understanding of the universal curvature jump at rough- 
ening transitions and facet edges and also provides an efficient way of locating 
the phase boundaries. 

PACS numbers: 68.35Rh, 82.65.Dp, 05.70.Fh 
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I. INTRODUCTION 



The theory of equihbrium crystal shapes (ECS) is well established Consider a 

macroscopic amount of solid in coexistence with and surrounded by its own fluid phase. 
The shape of the solid region is obtained by minimizing the total free energy of the solid- 
fluid interface subject to the fixed-volume constraint. This leads to the Wulff construction 
for ECS. Especially if one focuses on a particular direction, say the z direction, the crystal 
surface is defined by the height z{r) of the surface with respect to the position r = {x,y) in 
a reference plane. If the surface with slope m costs a free energy a{m) per unit base area, 
the ECS is given by H 

A^(r) = /(-Ar) , (1) 

where 2 A is the pressure difference between the two phases and /(h) is the Legendre 
transform of a{m): 

/(h) = min {o"(m) — h ■ m} . (2) 

{m} 

Here, h is the surface-tilting field conjugate to the surface slope. Equation (|l|) states that 
the surface free energy as a function of the surface-tilting field is itself the height of the 
surface from the base plane up to appropriate scaling of coordinates. 

Using this connection, the thermal evolution of the equilibrium shape of a face-centered- 
cubic (fee) crystal or a body-centered-cubic (bcc) crystal has been studied through the body- 
centered solid-on-solid (BCSOS) model which is equivalent to the six- vertex model The 
surface slope and the surface-tilting field in the BCSOS model correspond to the polarization 
and the electric field in the six-vertex model, respectively. The six-vertex model displays sev- 
eral ordered phases with ferroelectric or anti-ferroelectric order and a disordered phase . 
The disordered phase is a critical phase. Its scaling behavior is described by the Gaussian 
model, and is parameterized by the Gaussian coupling constant g or the stiffness constant 
K 0. At zero fields a Kosterlitz-Thouless (KT) type roughening transition takes place [§. 
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Below the roughening temperature Tr, the system is ordered with zero polarization when 
the field is small but it becomes rough with non-zero polarization beyond a critical value of 
the electric field. This transition is in the Pokrovsky-Talapov (PT) transition universality 
class [[7|,p!0|. The PT transition is characterized by the specific heat exponent a = 1/2 ||10|| , 
which implies that the free energy scales as / ~ |h — hd^^^. 

Equation (||) enables one to identify the ordered and the disordered phases to facet and 
rounded regions in ECS, respectively. Below Tr a facet appears surrounded by rounded 
vicinal surface. The rounded region in the ECS is rough in the sense that the height- 
difference correlation function behaves as 

{{z{r) - z{r')f) ~ In |r - r'l . 

The PT transition line corresponds to the facet edge and the crystal profile near the facet 
edge is given by z ~ {x± — x_Lq)'^/^, where x_l is a coordinate perpendicular to the facet edge 
and a;_L„ is the facet-edge position. 

A measurable quantity of physical importance is the surface curvature k, = 



^y"^, where the subscripts denote partial derivatives. From Eq. (|1|), the sur- 
face curvature is related to the Hessian of the free energy as 

A 



K = A VH [/(h = -Ar)] = . (3) 

'H[(T(m)J 



The Hessian of a function F(x) is defined by 

H[F(x)] = det 



F F 

^ x-i,xi ^ Xl,X2 



F F 



The second equality in Eq. follows from the identity H[/(h)]H[(T(m)] = 1. It has been 
predicted and measured experimentally that the curvature displays a universal jump 
at the roughening transition with discontinuity 

(A«:)kt = -77^ , (4) 
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where d is the distance between the crystal planes. It is also expected [0 that, for T < Tr, 
there is universal curvature jump at the facet edge with discontinuity 

(A«:)pT = . (5) 

These universal jumps are attributed to the universal nature of the roughening transition 
and the PT transition. Surface fluctuations without surface-tilting field are assumed to be 
described by the Gaussian model and the universal jump (Ak)kt is related to the jump 
of the stiffness constant at the KT transition [|ll|]. On the other hand, fluctuations of the 
vicinal surface near the facet edge are described by the one-dimensional free fermion model 



where the world lines of fermions are interpreted as step excitations in the surface [|T3[. It 
explains the universal jump (A/t)pT at the facet edge. 

In a recent paper on the six-vertex model, an exact relation has been found between 
the Hessian of the six-vertex model free energy and the Gaussian coupling constant g in 
the rough phase. The relation is given as H[/(h)] = (;^)^, where H is the Hessian of the 
six- vertex model free energy in units of ksT with respect to the dimensionless surface-tilting 
field. When one restores the dimensions, this relation becomes 

The Gaussian coupling constant determines the scaling exponents of various correlations 
and controls the finite-size-scaling (FSS) behaviors, i.e., a set of excitation energies AE of 
the transfer matrix, defined for a system with strip of width A^, satisfies the FSS form 

Re(A£) = ^^- + «^+Ar + ^j, (7) 

where C," is the imaginary part of the anisotropy factor, (A/", M) a non-negative set of integers, 
and (m, n) the level index, which takes integer values (see PJT^ for details). If one combines 
Eqs. dl) and (H), one gets 

2\cP 1 
hkbT g 



It was found that g = 1 at the KT-type roughening transition in the fee (110) surfaee and 
g = 2 at the PT-type facet edges If one uses them in Eq. the universal jumps in 
Eqs. (|) and (||) are obtained. Relation (H) or equivalently (|) is quite general in the sense 
that it does apply to the entire rough phase as well as at the phase transition points of the 
six-vertex model. A natural question that arises is whether such general relation is universal, 
in other words, whether it holds in other model systems too. This paper addresses to this 
question, and the results are affirmative. 

In Sec. H we introduce general models for crystal surfaces which display phase transitions 
between fiat and rough phases, and present renormalization group (RG) arguments that the 
Hessian of the free energy is the scale-invariant quantity. It yields the relation between the 
Hessian of the free energy and the Gaussian coupling constant in the rough phase. The result 
is given in Eq. (pTj). Combining it with the theory of the ECS, we obtain the relation between 
the surface curvature and the Gaussian coupling constant. To check the general theory, we 
present numerical results for the restricted solid-on-solid (RSOS) model in Sec. The 
Gaussian coupling constants, obtained in the FSS study from the transfer matrix spectra, 
is compared with the value obtained from the Hessian of the free energy. This confirms the 
validity of the universal relation. In Sec. ^V\, we discuss implications of the results and give 
brief summary. 

II. RENORMALIZATION GROUP THEORY 

Consider a solid-on-solid (SOS) type model for a two-dimensional crystal surface, where 
the surface is defined by height Zi at each site i on a substrate of size Li x L2 parallel to one 
of its crystal planes. The sites consist of projections of all lattice points on to the substrate 
and form a two-dimensional lattice. Figure 1 shows three examples of such substrates for 
sc (001), bcc (001), and sc (111) surfaces, respectively. In SOS type model the height is a 
single-valued function; there are no overhangs. The height at a given site can change by an 
integer multiple of the lattice constant 03 in the z direction. Due to the crystal structure 
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there may be p distinct classes of crystal planes parallel to the substrate, with inter-plane 
spacing d = a^/p. In that case substrate sites are separated into p sublattices and Zi takes 
the values {lid + integer x a^) if the site i belongs to the kth sublattice {li = 1,2, ... ,p). For 
examples shown in Fig. 1, p = 1, p = 2, and p = 3 in Fig. 1 (a), (b), and (c), respectively. 
Fee (110) surfaces also have p = 2. 

At low temperatures the surface will be in a flat phase. Steps which separate domains 
of fiat regions are the basic excitations and thermodynamic properties of the surface are 
described by a general Hamiltonian Ti which consists of the step-creation energy (?^s), the 
interaction energy between steps (Tii), and the surface-tilting energy (Tix) which controls 
the average slope of the surface. They are given by 

ns = JY.\^i-^j\' (9) 

Ht = -h{Az),L2 - h2{Az)2L, , (10) 

where denotes the pair of nearest neighbor sites, J is the step energy, (Az)i {{Az)2) 
is the total difference in height across the lattice in the 1 (2) direction, i.e., the difference 
in height between two sites i = {Li,y) and i = {0,y) {i = (x, L2) and i = (x, 0)) of 
the substrate, and h = {hi, /12) is the surface-tilting field. The explicit form of Tii is not 
important in the following analysis. Since long-wavelength fluctuations are dominant in 
the rough surface, we neglect height fluctuations inside the unit cell of the two-dimensional 
lattice and introduce a coarse-grained height zj which is the average height inside the unit 
cell containing sites i around the site j. (The coarse-graining scheme for the sc (111) surface 



is given in Ref. [|T5l.) Then {zj} takes integer multiples of d. 

The free energy as a function of the surface slope m = (mi, 7712) per unit base area is 
given by 

with the partition function Z = Y,'{zi} where the prime denotes the sum over 

all surface configurations satisfying the shifted boundary condition (SBC) {Az)i = iriiLi 



(i = 1, 2). m and h are related as h = VmO'(m)- Using the Poisson sum formula 

one can replace the discrete sum over z by integral over continuous field (j) with additional 
harmonic terms. The partition function in the continuum limit is then put in the form 

Z = f [V<P] e-'^«^ff['^Wl , (12) 

JSBC 

where the effective Hamiltonian is given by 

/9Hefr = | / dh \V<P\^ + pn' (13) 

with 

m' = -T.ynj rf'rcos (^^![!^^ ■ (14) 

The functional integral in Eq. (|12|) is taken over the field satisfying the SBC 

0(r + LiBi) = 0(r) + TTiiLi, (15) 

where ei (62) is the unit vector in the 1 (2) direction. Here we assume for simplicity that 
the substrate is a square lattice whose lattice constants in 1 and 2 directions are the same 
so that the stiffness constant K = 2(3J is a scalar. We will discuss later the more general 
case where the stiffness constant is a tensor. The sine-Gordon (SG) terms cos(27rn0/(i) 
account for the discreteness of heights. 

The effective Hamiltonian ?-^eff is the starting point of our RG arguments for the universal 
relation between the Hessian of the free energy and the Gaussian coupling constant. The 



RG theory for Ti^s with m = and in the absence of Tii is well established |T^-|T8[. At 
high temperatures where the surface is rough, all SG terms are irrelevant and hence V^'s 
renormalize to zero and K renormalizes to a fixed-point value K* . And at low temperatures 
the leading harmonics becomes relevant and the surface is fiat. At the roughening transition 
K* takes the universal value of |. The RG transformation for Ticfr in the presence of 
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Til and nonzero m can be performed similarly in the following way: (i) First introduce a 
field variable (f)o{r) = (f){r) — m • r, which satisfies the periodic boundary condition (PBC) 
(f)o{r + LiBi) = 0o(r). The partition function is then re-written as a functional integral over 
the field 0o as 

(ii) The field 0o is expended in a Fourier integral as 

J|p|<A (iTTj^ 

where A is the ultraviolet cutoff, and is separated into two parts 0o(r) and 0o(r) such that 
0o(r) = 0o(r) + 0o(r) and 0o (0o) only < |p| < A' (A' < |p| < A) components of the 
Fourier modes. Then the partition function is decomposed as 

JPBC 

j^jV<Pl] exp {-| j rf^rl V0'of - PU'Woij) + m ■ r + 0[;(r)]} . 

(iii) A partial integration over the fluctuations of 0o is performed and the remaining field 
00 is transformed back to 0'(r) = 0o(r) + m ■ r, which corresponds to the long- wave-length 
fluctuation part of 0(r). (iv) The RG transformation is completed by rescaling the momenta 
or the coordinate and the field as 

p— i>6p or r— i>r/6, 

0'(r) C0ne»(r/&) (16) 

with the scale factor h = A/A'. The scale factor ( for the field will be taken to be 1 to 
describe the Gaussian fixed point for the rough phase. Combining Eqs. (|15]) and (p!6D, one 
can see that 4>new satisfies the SBC 

(f>new{r + yBi) = 0„,e«,(r) + bnii (^y^ , (17) 

which implies that the slope is renormalized to 
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m' = 6m . (18) 

Under the RG transformation the free energy is transformed as 

a(P,m) = rV(P',m') + G , (19) 

where V denotes a set of model parameters, G is the analytic background depending on 
V and possibly on m, and V is the set of the renormalized model parameters. Let us 
focus on step (iii) where the partial integration over 0q is performed, which results in the 
renormalization of model parameters. However, the difference between the m = case is 
that the argument of Ti' is replaced by 0q — * {(p^ + m ■ r), which does not participate in 
the integral. Therefore K and the functional form of Ti' renormalize in the same way as 
at m = 0, and the renormalized values of V and G are only a function of V independent 
of m. As a consequence, one can readily see that the Hessian of the free energy is the 
scale-invariant quantity: 

H [a(m)] = H [a (m')] . (20) 

After successive applications of the RG transformation infinitely many times, the Hamil- 
tonian is renormalized to 

m* = lJ dh K'^IV^I' (21) 

with the renormalized stiffness constant K* provided the surface is in the rough phase. When 
the Hamiltonian is given by Eq. (^I]) , the slope dependent part of the free energy is easily 
isolated to be K*\m\'^ /{2j3) from a transformation 0(r) 0(r) — m ■ r so that the Hessian 
of cr(m) is simply given by {K* / (3Y. Thus from the scale-invariant property in Eq. (pOD, the 
Hessian of the original system is also given by 

H[a(m)] = [j^j . (22) 

In general, the stiffness constant in Eq. (|I3|) may be a tensor -ft'a,/? {a,/3 = 1,2). Following 
the same analysis, one can easily find that V and G in Eq. ( P^ ) do not couple to m either 
and Eq. (^) is replaced by 
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(23) 



where K^^ ^ is the fixed-point value of Ka^p. The Hessian of the free energy is also obtained 
from the scale-invariant property, which yields that 



H|a(m)] = . (24) 



The stiffness constant is not a good quantity since it depends on the scale of the field (j). 
So it is convenient to use the Gaussian coupling constant g which is defined as the coupling 
constant of the Gaussian model with the Hamiltonian 

PnG = ^ldh (25) 



where the periodicity of the field (p is set to 2n |[T9|. The periodicity of the field (f) is a^- So 
it is converted to 27r by rescaling ip = 2T[(f)/a^. After a rotation and rescale of coordinates, 
the Hamiltonian ( pB]) is transformed to the form of Eq. (^) with the Gaussian coupling 
constant given by 

g = 2^^det (K;^) (^y. (26) 

Using Eq. ( pG]) in Eq. ( pi] ) and H[/(h)] = l/H[(T(m)], one obtains that the Hessian of the 
free energy is given by the Gaussian coupling constant as 

' ^2 



H[/(h)] 



1 



(27) 



ksT J 2Tig 

The exact result of the six-vertex model in Eq. (^ is recovered since 03 = 2d in that case. 

Combining Eqs. (^ and (P7|), one finally obtains the universal relation between the 
surface curvature of the ECS and the Gaussian coupling constant in the entire rough phase 

^^^^ ^''^ 



where 03 = pd is used. Equations (|27| ) and (|28| ) are the main results of this paper. For fee 
(110) surfaces there are two equivalent crystal planes {p = 2). Inserting p = 2 in Eqs. ( P7[ ) 
and (|28|) , one reproduces the exact results of Eqs. (||) and (||). 



10 



III. NUMERICAL STUDIES OF THE RSOS MODEL 



In the previous section, we presented RG arguments for the relation between the Gaussian 
couphng constant and the Hessian of the free energy. It is obtained from the observation 
that the Hessian of the free energy is a scale-invariant quantity. In this section we test the 
validity of Eq. (P7| ) in the RSOS model on a square lattice (denoted by C) by comparing 
the Gaussian coupling constant obtained from the FSS amplitudes of the transfer matrix 
spectra and the value obtained from the Hessian of the free energy, using Eq. (P^). 

The RSOS model describes the surface of sc crystals viewed from the [001] direction. 
The Hamiltonian for the RSOS model with the surface-tilting field h = [hi, is given by 

^Rsos = K^5 {\zi - Zj\-l)- hi Y^{zi+e^ - Zi) - h2 ^^(^1+^2 - ^i)-, (29) 

where Zi is the integer- valued height variable at site i va. K is the step energy, 
denotes the pair of nearest-neighbor sites, and 61(62) is the unit vector in 1 (2) direction. 
The height differences between nearest-neighbor sites are restricted to and ±1. (In this 
section, length and energy are measured in units of lattice constant and /c^T, respectively. 
So all quantities are dimensionless.) The RSOS model with h = displays a roughening 



transition at = ~ 0.633 and g = 1/4 at the roughening transition [20]. There is one 
equivalent crystal plane parallel to the (001) surface. This means that the RSOS model 
represents a p = 1 case among the general cases discussed in Sec. |I|. 

Height configurations of the RSOS model can be mapped to arrow configurations on 
bonds of the dual lattice denoted hj Cd. If there is no step across a bond in Cd, no arrow 
is assigned to the bond. And if there is a step, an arrow is assigned in such a way that 
the height at the right-hand side of the arrow is higher than the other side by 1. Since 
there are no dislocations, the number of inward and outward arrows at each vertex should 
be equal (the so-called ice rule). There are nineteen vertex configurations satisfying the ice 
rule. So the RSOS model is equivalent to the 19-vertex model. The vertical (horizontal) 
slope corresponds to a net imbalance between left-right (up-down) arrows. 
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A row-to-row transfer matrix T is easily constructed. The partition function on a lattice 
of size N X M can be written as Z = Tr T^^. If one uses the PBC for the arrow variables 
in the N direction, the net number Q of up arrows in each row of vertical bonds is the same 
in all rows due to the ice rule. So the transfer matrix is separated into blocks of the form 

Q 

where Tq operates on the Qth sector, defined by the set of arrow configurations with net 
number Q of up arrows {Q = ~N, —N + 1, . . . , A^). The largest eigenvalue of Tq will be 
denoted by Ag = exp[—EQ{N)] and Eq{N) will be called the ground-state energy in the 
sector Q. In the limit N,M oo, the free energy e(mi, as a function of the horizontal 
slope mi = Q/N and the vertical surface-tilting field /i2 is given by 

e(mi,/i2) = hm . 

It is related to the free energy /(h) through the Legendre transform 

/(h) = min {e{mi, — hirrii} 

— l<mi<l 

and equilibrium values of mi and hi are related by hi = de{mi, h2)/dmi. 

It is well known that for a rough surface with average horizontal slope mi the transfer 
matrix spectra follow the FSS form 

EQiN) = Ne{mi,h2)-^, (30) 
Eq^^N) = Eq± nhi{mi, h2) + (31) 

where Q = rriiN, c = 1 is the central charge for the rough phase, (" is the imaginary part of 
the anisotropy factor, and g is the Gaussian coupling constant [^]. To obtain an estimate for 
g, one has to know the value of (" . It can be obtained from Eq. ( |30D using the two ground 
state energies Eq{N) and Eq/{N') for two values of strip width A^ and A^' chosen to satisfy 
the conditions Q = rriiN and Q' = rriiN'. Combining Eqs. ( PU] ) and (0), one can obtain 
the following estimate gFss{N) for g 
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9Fss{ ) ^2ArAr' 1^ N'Eq{N) - NEq^N') )' 
On the other hand if the relation (|27|) holds, it can be obtained from the relation 



^H(iV) = , ^ (33) 

2vrv/H[/(h)] 

as well, where H is the Hessian of the free energy in dimensionless form. The Hessian of / 
is directly obtained from the partial derivatives of e: 

H[/(h)] = . 

Cmi ,mi 

The partial derivatives are evaluated numerically as 

^^=NiEQ^, + EQ_,-2EQ), 

d^e _ {EQ{h2 + dh) + Egjh - dh) - 2EQ{h2)) 
dhl N{6h2f 

where we choose 5h2 = 0.001. This procedure gives estimates gn{N) for g. We use the 
subscripts in g to show how they are obtained. 

Estimates for g obtained in these two ways are shown in Fig. 2, where a data point 
represents a pair of values (gu^N) , gpssi^)) . Figure 2 (a) shows the results for mi = 
h2 = 0.0, and K = 0.2,0.4, and 0.6. For mi = 0, ^fss(A^) is obtained from Eq. (H) by 
choosing Q = Q' = and A^' = — 1. For each K, data shown are for the strip widths 

= 4, 5, . . . , 10 from left to right. A similar plot is shown in Fig. 2 (b) for mi = 1/2, 
K = 0.4, and /i2 = 0.1,0.15,0.2,0.25, and 0.3. For these cases, A^' in Eq. (H) is chosen 
to be A^ — 2 and the strip widths A^ are 6,8, and 10 from left to right. In all cases, they 
converge to the same values, i.e., the data points approach the line gu = gpss denoted by a 
broken line as A^ increases. The inset in Fig. 2 (a) shows the estimates for mi = h2 = 0.0, 
and K = ln[(v^ + l)/2], where the exact value of g is known to be 1/5 from the self-dual 
property of the RSOS model [^. Both quantities converge excellently to the exact value. 
The fact that g^i and gpss converge to the same value implies that the relation ( P7D holds 
in the RSOS model. Furthermore, as can be seen in Fig. 2, g-^ shows better convergence 
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than gFss- In addition to the better convergence property, gu provides a more convenient 
way of estimating g than the standard FSS study of the transfer matrix spectra. To obtain 
the estimate for one should have two ground-state energies at different strip widths with 
the same value of mi = Q/N. But it is difficult to find a set of integer values of and Q 
which gives the same value of mi. On the other hand, to obtain g^, one needs to evaluate 
the largest eigenvalues for a single value of A^. So the relation (^71) presents an efficient and 
convenient method to study the scaling behavior of the rough phase. 



IV. DISCUSSIONS AND SUMMARY 

In this paper, we introduce general models for crystal surfaces and derive, using the 
RG arguments, the relation (^) between the Gaussian coupling constant which determines 
the strength of critical fluctuations of the rough surface and the Hessian of the free energy. 
Combined with the theory of the ECS, it relates the surface curvature of the rounded region 
to the Gaussian coupling constant. In particular when applied to the phase transition 
point, it explains the universal curvature jump (Ak)kt at the roughening transition. The 
roughening transition takes place when the SG term in Eq. (|13]) becomes relevant. The RG 
calculations for that Hamiltonian shows that it has a scaling dimension Xp = p'^/{2g) [[TBI , 
which becomes 2 at the roughening transition, i.e., g = g^T = So one obtains the 

universal curvature jump 

2 Arf2 p2 2 \£ 
A/s: KT = - TT^J^ = -TT^- 34 

As examples, the RSOS model {p = 1), the BCSOS model {p = 2), and the triangular- 
Ising solid-on-solid (TISOS) model {p = 3) |21[] have the Gaussian coupling constants g = 
1/4, 1, and, 9/4, respectively, at the roughening transition points [^,^,0 and hence in each 
the universal curvature jump by the same amount. 

Below the roughening transition, there appears a facet which is separated from the 
rounded regions by the PT transition line. Near the PT transition systems become extremely 
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anisotropic and no conventional RG theory has been developed for the value of g at the 
transition points. Instead, the surface near the PT transition is studied using a random walk 
or free-fermion model [0, which predicts the universal curvature jump (A/t)pT in Eq. (^. 
Combining it with Eq. (pHD, one can see that the Gaussian coupling constant should be 
(yfpT = 2(yfKT at the PT transition points. Using that property, the PT transition point can 
be located accurately. In the RSOS model case, it is expected that = 1/2 at the sc (001) 
facet boundary. In Fig. 3 we present gFss{N) and gn{N) for a surface whose horizontal 
slope is fixed to (mi = 0) for several values of h2 below the roughening temperature. 
The surface remains fiat below a critical value of the surface-tilting field he. Above he the 
surface becomes tilted rough. The critical value of /i2 can be accurately determined from the 
condition that (7 = 1/2 at the transition. The insets of Figs. 3(a) and (b) show the estimates 
hc{N) for he obtained from the condition gpssi^) = 1/2 and gniN) = 1/2, respectively. 
Like the Gaussian coupling constant, hc{N) obtained from g-^ have less ESS corrections than 
those from gpss- The critical value h^ at K = 1.0 is estimated as 0.30 ± 0.01 (marked by 
arrows in the insets of Eig. 3) using polynomial fitting in j^, which is consistent with a value 
obtained from an alternative way |2^ . 

In summary, we derived the universal relation between the Hessian of the free energy 
and the Gaussian coupling constant in the rough phase of general surface model using 
RG arguments. It relates the surface curvature at the rounded region of the ECS to the 
universal quantity. Especially if it is applied to the phase transition points, it gives a clear 
understanding of the universal curvature jumps. The validity of the relation is checked in 
the RSOS model numerically. Erom the numerical results, it was found that the values 
of the Hessian have less finite-size corrections than the scaling dimensions obtained from 
the standard ESS theory. So, in practical points of view, this fact provides a better way 
in studying the scaling behaviors of the rough phase and the phase transitions in crystal 
surfaces. 
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Figure Captions 

Figure 1 Shown are the sc crystal on a substrate parallel to the (001) plane (a), the bcc 
crystal on a substrate parallel to the (001) plane (b), and the sc crystal on a substrate parallel 
to the (111) plane (c). The projected lattice points on to the substrates form 2-dimensional 
lattices and their sublattice structure is shown. 

Figure 2 The Gaussian coupling constants obtained from two different methods at mi = 
h2 = 0, and K = 0.2, 0.4, and 0.6 (a) and at mi = 1/2, K = 0.4, and /i2 = 0.1, 0.15, 0.2, 0.25, 
and 0.3 (b) are compared. The data are obtained from numerical diagonalizations of the 
transfer matrix for strip width = 4, 5, . . . , 10 in (a) and = 6, 8 and 10 in (b). The inset 
in (a) shows gu and gpss a-t mi — h2 — 0, and K — ln[(-\/5 + l)/2] where the exact value of 
g is known to be 1/5, whose location is indicated by the arrow. They show the converging 
behaviors to the exact value. The lines are guides to eyes. 

Figure 3 The Gaussian couphng constants gpss (a) and gn (b) are shown at mi — 0.0 and 
K — 1.0 > Kc ior several values of /i2- The insets show the estimates for the critical values 
of /i2 which are obtained by solving g = 1/2 numerically. The extrapolated values for he are 
marked by arrows. The lines are guides to eyes. 
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